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We investigate the chiral phase transition of a QCD-like theory, based the shape change 
of the effective potential near the critical point. The potential is constructed with the 
auxiliary field method, and a source term coupled to the field is introduced in order to 
compute the potential shape numerically. We also generalize the potential so as to have 
two independent order parameters, the quark scalar density and the number density. 
We find a tri-critical point locating at (T, fi) = (97, 203) MeV, and visualize it as the 
merging point of three potential minima. 



1. Introduction and Summary 

The phase structure of quantum chromodynamics (QCD) at finite temperature T 
and quark chemical potential \i is actively investigated in the context of the col- 
lider experiments using ultra-relativistic heavy ions, which create the highly-excited 
QCD matter.^ The chiral symmetry of QCD, which is spontaneously broken in the 
vacuum, will be restored at sufficiently high temperature and/or quark chemical 
potential. We investigate here the chiral phase transition of a QCD-like theory, fo- 
cusing on the shape change of the effective potential near the critical point. The 
QCD-like theory is the renormalization-group ( RG) im proved ladder approximation 
for the Schwinger-Dyson (SD) equation of QCDpEEl This theory describes the dy- 
namical chiral symmetry breaking, retaining the correct high energy behavior of 
the quark propagator. Using this theory and the low-density expansion, the pion- 
nucleon sigma term and the quark condensate at finite density have been calculated 
successfully.^ 

As for the effective potential, the Cornwall- Jackiw-Tomboulis (CJT) potential 
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functional has been used in several works for the study of dynamical chiral sym- 
metry breaking. However, the interpretation of the CJT potential away from the 
extremum is not obvious. Thus, in order to investigate the nature of the chiral 
phase transition through the shape change of the effective potential, we construct 
the Landau potential of the QCD-like theory^ by using the auxiliary field method, 
in which the bilocal external field is introduced. We generalize the potential so as 
to have two independent order parameters, the quark condensate and the quark 
number density at fixed T and ji. By plotting the 2D contour map in the order 
parameter plane, the shape change of the potential near the critical point can be 
elucidated. We find a tri-critical point is located at (T, [i) = (97, 203) MeV, and 
visualize it as the merging point of three potential minima. 

2. Basic Ingredients of a QCD-Like Theory 
2.1. In the vacuum 

In the QCD-like theory, first, we neglect the self-interaction between the gluons 
and integrate out the gluon fields, leaving the non-local interactions between the 
quark. The non-Abelian nature is taken into account in the running of the cou- 
pling constant. The generating functional Z = J VipT>'t/jT>Aexp(i J c1 4 x£qcd) is 
approximated by Z~f DipDip exp(i J d xC) where 



Jp 

~\ I p j a ( p ~ !) + 1) Ka " nS (p> k ^{ k + f) M k - 1) -w 

Here, ip(p) is the quark field in the momentum space and J = Jd 4 p/(27r) 4 . The 
indices, a, /?, • ■ • correspond to the Dirac structure. We defined the kernel K as 

K a ^ s {p, k) = g 2 ( lll T a ) Sa ( ly T a f-<iD^{p - k) , (2) 

where T a is the color su(N c ) generator and iD^ v (p) is the gluon propagator: 

iD^(p) = 9 - [ P /P . (3) 

Hereafter, we employ the Landau gauge, a = 0. The non-Abelian nature of the gluon 
interaction is treated in this model as the one-loop running of g: g{pa.ax(j>%, k%)) 
with pe the momentum in Euclidian space. The divergence of ~g(p%) appearing at 
Pe = Aqcd is removed by introducing an infrared cutoff parameter pie as^ 

2 1 

f {pl) = a MCpWfVA^cd) (4) 

with a = 1/(8tt 2 ) • (11N C - 2iV»/3. Here, N c (= 3) and N f (= 2) are the numbers 
of colors and flavors, respectively. 
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We next introduce the following bilocal auxiliary field for non-perturbative anal- 
ysis: 

Xa/J(p, q) = J K a ^<\ Pl k)^ (k + |) ^- (k - |) . (5) 

This makes the action bilinear in the quark fields ip and ip. Then, integrating out 
ip and ip, we end up with Z = J Pxexp? T[\] where the classical action reads 

r[x] = \t tT[x(p,-q)K- 1 (p,k)x(k,p)}- l TTLn(^ i (q)(2n) i - X (p,-q)).(Q) 

^ J pkq 

In the mean-field approximation the SD equation is obtained as the extremum 
condition for this classical action for x- The existence of a non-trivial solution for 
X indicates the dynamical breaking of the chiral symmetry. In the vacuum this 
solution is expressed as (Xap) — E Q/ 3(p)<5 4 (g)(27r) 4 with the mass function £, and 
the effective potential V[S] = — T[x]/ J d 4 x becomes 

V\Z] =-\f tr[E(p)lf- 1 ( P ,fc)S(fc)]+i f tr!n(i/-S(p)). (7) 

^ Jpk J p 

Then the SD equation, SV/5T, = 0, in the improved ladder approximation is written 

Although the E aj g(p) has the general form T,(p 2 )5 a i3 + £t,(p 2 y Q(3 in the vacuum, 
we can set Y> v (p 2 ) = in the Landau gauge (a — 0). After carrying out the Wick 
rotation and the angle integration, the SD equation becomes 

s <»-> = ^ f ^(■^(p|.4)) mM( 4 -4) k% 1%, ■ m 

where C 2 {N C ) = T a T a = (A 2 - 1)/2JV C . 



2.2. Order parameters 

The quark condensate with the four momentum cutoff A is defined as 

This bare value at the scale A is converted into the value at the lower energy scale 
\i (e.g., 1 GeV) via the renormalization group equation 

l/[4((12C 2 (iV c ))- :l (llJVc-2iV / )/3)] 



(ipip) u = W)a ( hrrx ) • ( n ) 



,5 2 (M), 

Next, the pion decay constant f v is estimated in terms of the mass function 
E(p 2 ) by utilizing the Pagels-Stokar formula: 

h ^ Jo Pe (pI + mvl)? V Pe) 2 dp% ) ■ (12) 

We fix the value of Aqcd here so as to reproduce the empirical value of f n . 
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2.3. At finite temperature 

We use the imaginary time formalism to extend the QCD-like theory to the case 
with finite temperature and chemical potential, making the following replacement: 

/ f{Po,p) > T jr I ^Pf(iuj n + u,p) , (13) 

J P n=-oo J ^ A7r > 

where u> n — (2n + l)nT (n G Z) is the Matsubara frequency for the fermion. 
The mass function ~E a p(uj n , p) at finite T and /i decomposes into E(w n , [p|)<5 a /3 + 
E s (w„, \p\)cJ n (lo)a/3 + S„(w„, |p|V(7i)a/3- We assume here E s = E„ = for sim- 
plicity for our purpose of demonstrating the usefulness of the effective potential in 
the QCD-like theory. Further, we use a covariant-like ansatJ^for the mass function: 

%,P)^S(P 2 ), (14) 

where the frequency and the momentum appear in the combination p 2 = uj 2 + \p\ 2 
in E. 



3. Landau Potential 

Here we explain a method that we use to construct the functional form of the ef- 
fective potential away from the extremum. A standard way to assess the potential 
form is to apply an external source which is coupled to the field linearly. It is rec- 
ognized, however, that using this approach we cannot study non-convex potentials, 
which is an important feature near the phase transition point. We thus apply an 
external source field J(p, k) which is coupled to the square of the self-energy as 

V\S, J] = V[E] + \ [ tr[E(p) J(p, fc)E(fc)], (15) 

z Jpk 

where V[E] has been defined in J7J. By imposing the extremum condition for V[E, J] 
with respect to E(p), we derive the SD equation with the source field J and have a 
solution Ej(p). The effective potential V for the configuration Ej(p) is written as 

V\Ej] = V\Ej, J]- \ ( tr[Ej(p) J(p, k)Xj(k)} . (16) 

1 Jpk 

Near the critical points we compute T^[Ej] for a family of Ej obtained through a 
particular set of the function J(p, k). Here, among infinite possibilities, we choose 
one natural choice, J(p, k) — —cK~ l (p, k), with a parameter c, and study the shape 
of the effective potential along this variation. 

In the numerical evaluation of the effective potential, we split V[E] as V[E] = 
(V[E] - V[0]) + V[0]. Here, the difference between the free energies V[Y,j] - V[0] 
is finite and can be evaluated numerically. The last term V[0] =i J tr In is 
the potential of a free massless quark gas and is divergent due to the vacuum 
fluctuations. In our regularization to remove this divergence, we replace the last 
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term V[Q] with the pressure Pf ree of a free massless quark gas: — V[0] — > -Pftoc = 
N c N f T A [7tt 2 /180 + (1/6) • (p/T) 2 + (1/12tt 2 ) • (p/T) 4 . 

In order to introduce another order parameter, the quark number density p, 
we first replace the variable p with the quark number density p through the usual 
Legendre transformation, 

F((^},p,T) = V((^),p,T)+pp (17) 

with p = —d V({if>tp), p, T)/d p. Then, by adding a coupling energy with an external 
potential v to Eq. (fT?]) . we define a new effective potential as 



V((W),p, v, T) = F((W),P, T) - pv, (18) 

which may be interpreted as a Landau potential!^ Note that p and v are indepen- 
dent variables here. When the condition dV /dp — pt,— v = is imposed, this new po- 
tential V coincides with the original potential: V(('ipij}}, p(p); p,, T) = V((i/jip); p, T). 



4. Numerical results 

We plot V[Y, C ] - V[Q] as a function of |(V>V)I/Aqcd in Fi S- !( a ) at A* = °> 100 and 
300 MeV for several values of T and in Fig. 1(b) at T = 80 and 100 MeV with 
several values of p. As seen in Fig. 1(a), a second-order transition occurs in the 




Fig. 1. The effective potential, V[S C ] — V[0], with (a) fixed fi and (b) fixed T as a function of 
\(j[np)\. (c) The contour map of the Landau potential in the p-fyrp) plane with various chemical 
potentials for T = 100 MeV (left) and T = 80 MeV (right). 
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range T — 100 - 120 MeV in the cases with p = and p = 100 MeV. At p, = 300 
MeV, however, a first-order transition occurs between T = 75 MeV and 80 MeV. 
Similarly, in Fig. 1(b), a first-order transition is seen to occur when the chemical 
potential is changed between p — 280 MeV and 290 MeV at low temperature 
(T = 80 MeV), while the transition is second order at fixed T — 100 MeV. 

In Fig. 1(c), we plot the contour map of the Landau potential given in (| 18|) 
at several values of p with fixed T, namely T = 100 MeV (left) and T = 80 
MeV (right). The vertical axis represents the quark condensate {-ipi/j} / Aq CB , and 
the horizontal axis represents the quark number density in units of the nuclear 
saturation density po- It is shown that a second-order transition occurs in the left 
panels of Fig. 1(c). At /z = 0, the two minima are positioned symmetrically in the 
p-{i[)ip) plane at p — with finite values of the quark condensate, (ipip)- As the 
chemical potential increases, the two minima approach each other for finite p, and 
then fuse continuously to form a single minimum at (i/)ip) = an d finite p. At lower 
T (in the right panels of Fig. 1(c)), a new local minimum appears at (iptp) = as p 
increases, and three local minima exist in a certain range of p. These three minima 
become energetically degenerated at the critical point, and we observe a first-order 
chiral phase transition. Note that a tri-critical point is located at (T t , pt) — (97, 203) 
MeV, where these three minima merge. 
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